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Abstract. We propose an approximation scheme to describe the dynamics of the 
spin-boson model when the spectral density of the environment shows a peak at 
a characteristic frequency f2 which can be very close (or even equal) to the spin 
Zeeman frequency A. Mapping the problem onto a two-state system (TSS) coupled 
to a harmonic oscillator (HO) with frequency wo we show that the representation 
of displaced HO states provides an appropriate basis to truncate the Hilbert space 
of the TSS-HO system and therefore a better picture of the system dynamics. We 
derive an effective Hamiltonian for the TSS-HO system, and show it furnishes a 
very good approximation for the system dynamics even when its two subsystems are 
moderately coupled. Finally, assuming the regime of weak HO-bath coupling and low 
temperatures, we are able to analytically evaluate the dissipative TSS dynamics. 
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1. Introduction 

As new experimental physics techniques allow us to reach minute length scales, and 
milestones in control and precision of the performed experiments are achieved, the 
study of open quantum systems has become one of the most prominent areas of modern 
physics. Moreover, both fundamental and practical aspects of the systems under current 
investigation help to foment a very broad interest in the area. On the fundamental side, 
the possibility of having a better comprehension of how classical physics emerges from 
a quantum system or the realization of very peculiar quantum mechanical states have 
attracted the interest of many researchers towards the area. From the practical point of 
view, understanding the main mechanisms leading to the loss of quantum coherence in 
such systems constitutes a key feature, and one of the major challenges for the physical 
implementation of, for example, a quantum computer. 

A very successful open quantum system which has been investigated to date is a 
two-state system (TSS) coupled to a dissipative environment. Despite its simplicity, 
the TSS dissipative dynamics is the paradigm of a wide variety of physical systems. 
Superconducting devices containing Josephson junctions pQ, few-electron semiconductor 
quantum dots [2] and two-level atoms in optical cavities [3] are just few examples of 
systems whose low-energy level dynamics can be, in general, captured by a dissipative 
TSS. 

The most studied model of a dissipative TSS is the well-known spin-boson model 
[I], where the TSS is linearly coupled to the coordinates of a bath of non-interacting 
oscillators. The dissipative effects of the TSS phenomenological environment are 
determined by the spectral density of the bath of oscillators [S], which, in general, 
is assumed to have a power-law behavior at low frequencies. Since it seems that no 
general analytical solution can be obtained for such a model, several approaches have 
been proposed to determine the time evolution of the dissipative TSS . Those approaches 
basically belong to two distinct approximation schemes: one considers a weak TSS- 
bath coupling in the low-temperature regime, and the other consists in performing 
an expansion in the tunnelling amplitude A of the TSS states. The first approach is 
commonly implemented by either using path-integral methods [H] or the Born-Markov 
approximation (also known as Bloch-Redfield formalism) [3, El E] , where only the lowest 
order terms in the TSS-bath coupling are taken into account for the TSS dissipative 
dynamics. The other scheme has been shown to be a fair approximation specially in the 
regimes of strong damping and/or high temperatures, and also employs path-integral 
methods within the so-called noninteracting blip approximation (NIBA) as its main 
technique jH E] . 

An emerging problem in the area of dissipative TSSs is the one in which the 
bath effective spectral density J e ff(u) presents a pronounced peak (resonance) at a 
characteristic frequency Q. A typical example of such a case happens when the energy 
scale of the device used to detect the state of the TSS is comparable to that of its 
own regime of operation. Since the device, or "quantum detector", also suffers from 
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the dissipative effects of the environment, the TSS-detector resonance constitutes an 
efficient channel for decoherence processes to take place. For example, as the energy 
scale of superconducting qubits reaches the order of several GHz, this mechanism shows 
up on these systems due to their coupling to the read-out dc-SQUIDs [TOl [IT]. In 
addition to the TSS-detector case, the presence of a structured bath with a sharp peak 
has also been discussed in other pertinent contexts, e. g., electron transfer in biological 
and chemical systems [12J and semiconductor quantum dots [13], to name just a few 
relevant problems of this growing interest subject. 

Indeed, it has been shown [TJ1 [151 E] that the structure of the environment 
becomes an essential feature of the dissipative TSS dynamics when its frequency is 
comparable to the bath resonance frequency fl Not only the decay rates cannot be 
correctly determined by a perturbative approach in the TSS-bath coupling, but such an 
approximation also fails to account for the presence of additional resonances in the TSS 
dynamics. By using an exact mapping [12] of the spin-boson problem with a structured 
environment onto that of a TSS coupled to a single HO system with frequency Q, 
which itself interacts with a bath of oscillators having an Ohmic spectral density, several 
techniques have been applied to the investigation of the TSS dissipative dynamics. When 
the single HO is strongly damped and/or coupled to a bath at high temperatures, the 
NIBA has been employed to both TSS-HO weak [TT] and strong [TS] coupling cases. 
On the other hand, the numerically exact method of the quasi-adiabatic propagator 
path-integral (QUAPI) [HI [TT], as well as a three-level system (3LS) Redfield theory, 
have indicated the failure of the perturbative approaches in the HO-bath weak coupling 
and low temperatures regimes. 

In this work, we present an approximation scheme to describe the dynamics of the 
spin-boson problem with a structured bath having a pronounced resonance at frequency 
Q. Following previous works [121 HH [T3], we map this problem onto a TSS coupled to a 
single HO with frequency ujq. In order to obtain an appropriate effective Hamiltonian for 
the lowest-lying energy states of the TSS-HO system, we show that the representation 
of displaced HO states whose displacement depends on the TSS state, provides a better 
picture of the system dynamics, and consequently the appropriate basis to truncate the 
Hilbert space of the TSS-HO system. We show that the derived effective Hamiltonian 
furnishes a very good approximation for a vast range of the physical parameters of the 
system even when the TSS and the HO are moderately coupled close to the resonance. 
In addition, assuming the regime of weak HO-bath coupling and low temperatures, we 
are able to analytically evaluate the dissipative TSS spin dynamics. 

This paper is organized as follows. In section [2] we present the model for the system, 
and derive the effective Hamiltonian using the approximation scheme proposed above. 
By varying the physical parameters of the problem, we compare its predictions with: 
(1) an exact numerical calculation (for which we found to be sufficient to consider 
a N = 18 dimensional TSS-HO Hilbert space); and (2) the simple truncation of 
the TSS-HO Hilbert space. Section [3] contains the results obtained for the TSS- 
HO dissipative dynamics assuming the regime of HO-bath weak-coupling and low- 
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temperatures. Finally, in section |4j we present our concluding remarks. 
2. The model 

Our starting point for describing the TSS dissipative dynamics in the presence of a 
structured environment is the spin-boson Hamiltonian |3], 

H S b = + ^tyfik + ft°z ^2 Ck $k + t>h), (1) 

k k 

where <7j are the Pauli matrices and HA represents the tunnelling amplitude between 
the TSS states. 

The bath of oscillators, introduced by the canonical bosonic creation and 
annihilation operators b\ and b k , is assumed to have a spectral density J e //(co>) showing 
a Lorentzian peak at a characteristic frequency Q. 

It has been shown [T21 D3J that the proposed problem can be mapped onto a system 
comprised of a TSS coupled to a single HO of frequency uo, which itself interacts with 
a bath of oscillators having a spectral density J{uj) presenting a power-law frequency 
distribution. The system Hamiltonian in this formulation is given by 

H = — ^Ao" x . + hjjQa)a + ha z {ga) + g*a) 

+ ^tm k b\h k + h{o) + a)^2\ k (b{ + b k ), (2) 

k k 

where a) and a are the creation and annihilation operators of the single HO system, 
and g stands for the TSS-HO coupling constant. If g is a real function, the TSS-HO 
coupling has the form of a "a 2 -coordinate" interaction. On the other hand, if g is a purely 
imaginary function, the TSS-HO coupling becomes a "^-momentum" interaction. 

From the system Hamiltonian (J2|, one can observe that the dissipative effects 
of the environment on the TSS dynamics occur only indirectly through the TSS-HO 
interaction. Thus, in this formulation, it is clear that the TSS-HO coupling works 
as a dissipative channel for the TSS. Such a channel can be enhanced or suppressed 
depending on how strongly coupled the TSS and HO systems are. This depends not 
only on the strength of the coupling constant g, but also on how far from resonance these 
two subsystems are. Indeed, if we assume a "<j 2 -coordinate" coupling, i. e., g — g*, and 
an Ohmic spectral density J(u>), 

J(w) = | J2 X l 6 ^ - = (3) 

k 

one can show [12] that the corresponding exact mapping between the Hamiltonians in 
({!]) and ([2])[|] occurs when the single HO has its frequency given by the J e ff(u>) resonance 
frequency Q, i. e., u = Q. Consequently, the peak in the effective spectral density of 

| Assuming a Lorentzian shape J e ff{u>) = ~ = (n 2 ~^^+(2tthujQ.) 2 ^ or ^ e spectral 

density of the system Hamiltonian (flj) . 
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the bath "seen" by the TSS, J e ff(u), can be interpreted as a manifestation of a definite 
resonance of one of the environment degrees of freedom. 

Recently, Westfahl Jr. and collaborators [13] have demonstrated that a u a z - 
momentum" TSS-HO coupling, i. e., g = —g*, with a power-law frequency dependence 
for J{oj) (J(u) oc u s ), leads to a significant change of the effective bath spectral density 
J e ff(u) of the system Hamiltonian ([T|. Although J e ff(uj) still shows a pronounced peak, 
featuring an approximated Lorentzian shape, its characteristic frequency Q does not 
match the single HO frequency u , i. e. fl ^ uo. Indeed, investigating the dissipative 
effects, due to the spin-orbit mechanism, on the electronic spin dynamics trapped in 
quantum dots, they have found that the characteristic resonance frequency Q predicted 
for the effective bath spectral density J e ff{uj) is, in general, shifted to a value much 
lower than that of the single HO, Q <C ujq. In addition, they have shown that such 
a resonance occurs in a frequency regime which can be, in principle, of the order of 
the Zeeman frequency of the spin, raising concerns about the appropriate approach to 
correctly describe the dissipative spin dynamics. 

Thus, since the system Hamiltonian ^ is capable of describing the dissipative 
dynamics of a TSS coupled to a variety of structured environments, from here onwards, 
we shall focus on describing its features and properties. Unfortunately, despite its 
simplicity, the Hamiltonian ^ cannot be analytically diagonalized. However, since our 
final goal is to determine the TSS dissipative dynamics in the regime of low temperatures, 
HA ^> ksT, we can concentrate on determining its low-energy level spectrum. In 
addition, we are going to assume that the HO dynamics is subject to weak dissipation, 
in such a way that its states are weakly perturbed by the HO-bath coupling. 

Under these considerations, in order to find a low- dimensional effective Hamiltonian 
for the TSS-HO system, the natural procedure is to perform a truncation of the HO 
Hilbert space, thus reducing the TSS-HO Hilbert space dimensionality. The simplest 
method of doing that consists in considering only the ground, |0), and first excited, |1), 
states of the HO, leading to a four- dimensional Hilbert space of the composite TSS-HO 
system. 

However, as depicted in Fig. [TJ that approximation fails as soon as the condition 
of a TSS-HO weak coupling, g <^ A,u , is not satisfied. Comparing with a numerical 
calculation performed by taking into account a N = 18 dimensional Hilbert space, 
we can see the system Bohr frequencies ujq\ and u m such that ftcuoi = E% — E and 
hw Q2 = E 2 — E , where E , E\ and E 2 respectively correspond to the eingenenergy of 
the ground, first and second excited states, can deviate circa 10% from the correct 
values, even for moderate TSS-HO couplings. Since we expect a sharp resonance for the 
effective spectral density J e ff(u), those deviations can easily lead to much larger errors 
when computing the TSS relaxation rates. Indeed, Thorwart and collaborators [13] 
report deviations for the relaxation rates of the order of ~ 20% for the weak coupling 
case g = 0.07, when performing the simple truncation of the HO Hilbert space. 

The failure of the simple truncation of the HO Hilbert space for cases other than the 
TSS-HO weak coupling, g C A, O , was already expected, since it does not correctly take 




Figure 1. (a) The four lowest levels of the TSS-HO system, as a function of the 
tunnelling frequency A, for the moderate coupling case g — 0.3S1, obtained through: 
numerical calculation considering a N — 18 dimensional Hilbert space (black solid 
curve) ; analytical solution after performing a simple truncation of the HO Hilbert space 
(blue dot-dashed curve); and analytical solution of the effective Hamiltonian derived, 
Eq. (J8j, in the HO displaced state representation (DSR) (red dashed curve). (The 
ground state is taken having E — 0). (b) and (c) present the Woi and W02 frequency 
deviations, AujQi = u>q^ ,DSR — u!q^~ 18 \ for the simple truncation approximation and 
the effective Hamiltonian Eq. |8j, as a function of the TSS-HO coupling g (for the 
resonance case A = f2), and the tunneling amplitude A (assuming g — 0.4S7) . 



into account the weight of the higher energy HO states in the spectral decomposition of 
the lowest energy eigenstates of the Hamiltoninan (J5]). As one can observe from the latter 
equation, the TSS-HO interaction term has the net effect of inducing displacements of 
the HO in its phase space {q,p}, hence leading to the mixing of the HO eigenstates of 
Eq.([2]). So, the stronger the TSS-HO coupling is, the more important the HO higher 
energy states become. 

Therefore, in order to obtain a better picture of the lowest energy states of Eq. 

we have to perform the truncation of the HO Hilbert space in a representation 
where such an effect is taken into account. In other words, it is the displaced harmonic 
oscillator states which ought to be the appropriate basis where one should perform the 
reduction of Hilbert space of the TSS-HO system. 

The change of representation can be done by introducing the unitary HO 
displacement operator D(at) = e aat - a * a ; which is very well-known in the context of 
coherent states [3]. There, among several other implications, the parameter a is a c- 
number, such that the transformation (a)aD(a) = a + a has the action of shifting 
the HO canonical position and momentum operators by amounts proportional to, Re(a) 
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and Im(ai), respectively. However, an important feature of the Hamiltonian ^ is the 
fact that the displacements of the HO states depend on the state of the TSS. Therefore, 
we are led to introduce the unitary transformation 

D(s,a z ) = e ^- s *^ % (4) 

as the conditional displacement operator of Eq.Q. Here, we consider the displacement 
parameter s as an ad hoc parameter of our model, to be properly chosen in order to 
obtain the best effective Hamiltonian possible. 

Transformation ^ is very similar to what is known in the literature as the polaron 
transformation [U [6]. Indeed, if we had followed, by analogy, the choice made for 
the polaron transformation, we would have obtained s = —g/uj^. Although it correctly 
diagonalizes the Hamiltonian ^ when A = 0, this choice is expected to fail, as discussed 
in [IT], when the tunnelling amplitude becomes appreciable, since from its analytical 
solution one obtains cj i — > oo, instead of u>oi — > uj , when A — > oo. 

Thus, in order to make the best choice for the ad hoc parameter s, we have to 
carefully observe which one would reflect the kind of dynamics we should expect for our 
TSS-HO system. Performing the change of representation, H = D^HD, and truncating 
the HO Hilbert space, one finds 

h ~ 

H T ss-ho = ~2^ x i 1 ~ 4 l s ! 2 c t c) + fkOQ&c 

+hc^ |<7z(<7 + uj s) + ifrysAj + h~c |<7js(<7* + uj s*) — ia y s*Aj , (5) 

where, in order to stress the fact that we are now working with a reduced Hilbert space, 
we have defined the operators c = |0)(1| and & = |1)(0|, in which |0) and |1) represent, 
respectively, the ground and first excited states of the displaced HO. A represents the 
renormalized tunnelling amplitude due to the TSS-HO coupling, which is given by 

A = Ae- 2|s|2 . (6) 

Without dissipation, the TSS-HO system should evolve in time in such a way there 
would be a persistent exchange of a quantum of energy between the two subsystems. 
Such a dynamical behaviour has a precise analogue in the problem of absorption 
(emission) of a single photon by (from) an atom placed in an optical cavity. There, 
the fundamental tool for the study of the system dynamics is the so-called Jaynes- 
Cummings (JC) Hamiltonian. Inspecting Eq. ([5]), one can see that the choice 

s - — g -rrr CD 

uj + |A| 

leads us to the same form as a JC Hamiltonian. It worth noting that one could use 
a variational calculation to determine the value of the parameter s for the polaron 
transformation Eq. (|4]). Indeed, by minimizing the free energy of the system using 
the Bogoliubov-Peierls bound, Silbey and Harris [18j could determine a transcendental 
equation for s. For the unbiased case and T = 0, that leads to a very similar renormalized 
tunnelling amplitude to that found in Eq. (fsj) (cf. Eq. (103) in the review |19j ) 
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Assuming A < C[§J the effective Hamiltonian in the HO displaced state 
representation (DSR) obtained using Eq.0 in Eq.Q reads 



Hdsr = ( l ~ 7 — ~^T7^ gtg l + ^o^c 

2 V (w +A) 2 J 



(wo + IAD 1 



j r (g&&™+g't&™), (8) 

Wo + |A| 

where we have introduced the ladder operators for the spin x— component, <r± = |(o" z T 
^!/]QIl The Hamiltonian ([8]) can be diagonalized analytically, yielding eigenenergies given 
by 



E = --\A\, E 3 = -\A\\l ^ — ^ + ^ , and (9) 




E12 = - < 2lo — ' =F W ( . — 2u + m r I . (10) 

The analysis of the obtained eigenenergies Eqs. ([9]) and (10), reveals that 
the Hamiltonian ([8]) provides the correct results for certain limits of the parameters of 
Eq. ([2]), namely: (1) as cu — ► 00, we obtain cjqi - ► |A|; (2) for |A| — > 00, we have 
uqi — > ujq; and for the resonant case A = ujq, with g <C A,^, we find ujqi pa ujq — g 
and ujq2 pa uj + g. Moreover, as one can see from Fig. [TJ the resulting Hamiltonian 
(Jsj) gives a very good approximation of our problem for a vast range of the physical 
parameters A and g. In fact, the largest deviations we have found considering the 
parameters illustrated in Fig. [I] were Acj 01 pa 0.07% and Alu 02 ~ 0.6%. It is worth 
noting that, when performing the simple truncation, this level of agreement can only be 
met for the ground and first excited states if we consider a N = 8 dimensional Hilbert 
space. Indeed, for a N = 6 dimensional Hilbert space we found: Au4f " 6) ~ 20Aoog SR 
and Auj { 2 =6) pa 3Au^ 2 SR , while for the case iV = 8 we obtained: Auj^ =8) pa 1.3Au^ SR 
and Au ( 02 =8) « 0MAu° sr . 

Finally, defining the states |±) as the eigenstates of a x , a x \±) = ±|±), the 
eigenstates of the Hamiltonian ([8| can be written in the form 



|e ) = |-,0), |e 3 ) = |+,I), \ ei ) = A\+,0)+B\-,l), and (11) 
\e 2 ) = -B\+,6)+A\-l), (12) 



where, A= (E + E 2 )U (E + E 2 ) 2 + A\g\ 2 A 2 h 2 /(uj + | A|) 2 and B 2 = 1 - A 2 



The eigenstates |ei 2 ), Eqs. (11 ) and (12), reveal the hybridization process occurring 



between the TSS and HO states due to their coupling. Indeed, if g = 0, i. e., the 
case of a decoupled TSS-HO system, the eigenstates of Eq. d8|) are simple tensor 



§ For the case A > 0, the last term of Eq. (JsJ> becomes (g$a^ + g*c a^), reflecting the fact that 
the state |+) has the lowest Zeeman energy. 

I The ladder a± operators obey the usual commutation relations: [a x , <Tjf **] = ±2<r±^ and 
[cr+\ <j^] — a x . Their action on the eigenstates of a x are given by: 0+ |— ) = |+) and eri 2: ' ) |+) = |— ), 
with (a^) 2 = {a^Y = 0. 
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product of the TSS and HO eigenstates (A = 1). Nevertheless, as one increases the 
coupling g, the system eigenstates become more hybridized. At the resonance condition, 
cuo = A, the eigenstates jei^) become nearly the maximally entangled TSS-HO states, 
with A ps ^=(1 + 0(g)). 

3. Dissipative dynamics 

The evaluation of the dissipative TSS dynamics can be described, for a weak HO-bath 
coupling and low temperatures, using the Redfield equations [6] 

Pnm(t) = —iUnmPnmit) — ^ Rnmklpkljt) , (13) 

k,l 

where the matrix elements p nm = (n\p\m) are taken using the eigenstates of Eq. p|, and 



the Redfield tensors are given by R nmkl = 6 lm £ r ^lk + 6 nk E r r LL ~ r !mL ~ r LL 
with 

POO 

r£i* = ^ 2 / dte-*-"**(W /tIm (0W/^(0)>, (14) 
Jo 

/■oo 

rLL = ^ 2 / dte^*(Hi M (0)Hi,r*(t)), (15) 
Jo 

where TCi(t) = e lHRt ^ h Hie~ tHRt ^ h is the interaction Hamiltonian in the interaction 
picture, and the brackets represent thermal averages over the bath degrees of freedom. 
Here, we denote Hr and Hi as the bath and HO-bath interaction Hamiltonians 
obtained from Eq. (|2]) under the transformation presented in the last section. The 
full Hamiltonian now reads 

H = H DSR + fokblh +(c ] + c- 2Re( ^, ^) £ hXk ^ + ( 16 ) 
k ^ ^o + |A| / fc 

where the second and third terms on the RHS of this equation are, respectively, H R and 



Hi. It is worth noting that the term proportional to Re(g) in Hi, Eq. (16), represents 
a correction naturally introduced by the scheme proposed, and leads to an effective 
interaction between the TSS and bath mediated by the TSS-HO coupling. Had we 
followed the simple truncation approximation, such a term would not be present in the 
system dynamics. We have seen that this term can represent a non- negligible correction 
to the system-bath interaction, since it can reach fractions of few percent of the latter, 
e. g., for g = 0.1, A = cu = 1, we find ||2Re(#)/(^o + I^DII/II^ + c|| ps 0.1. 

In order to be able to analyze the resonance case, u = A, as pointed out in 
Ref. [Hj, we have to go beyond the two-level approximation when using the Redfield 



equations (13). Indeed, since a two-level approximation is not capable of taking into 
account the hybridization of TSS-HO states, it fails to present the additional resonances 
of the TSS dynamics. Thus, considering a three-level system, one can determine the 
expectation value of a z , a z = Tr(p(t)a z ), as 

a z (t) ps 2Re [A{e \p\e x ) - B(e \p\e 2 )) . (17) 
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Performing the secular approximation [7J, we can neglect the terms R nm ki of p nm for 
which uj nm 7^ Ufa- This leads us to the following set of coupled differential equations 

An = ( e o|p|ei) « iwoiPoi ~ 7iPoi - 72P02, (18) 
P02 = ( e o|p|e 2 ) « ^02^02 - 73P01 - 74P02, (19) 

where 71 = #0101,72 = #0102,73 = -R0201, and 74 = #0202- Assuming that J(u) is a 
regular function in the complex plane, we can analytically evaluate the Redfield rates, 
obtaining 

Re T lmnk = (f^)J(M^ mu ^ /2) (20) 
Im V\^ nk = --(h lm h mk )P [ du 7^ 2 \0J - uj nk coth(^J^) ) (21) 

71 J UJ - LO nk \ I J 

where h lm = + c - ^^d: z )\m) . In addition, using the relation Y { ~^ ml = (rf^)*, 

we can determine the Redfield tensors R nm ki in terms of the rates r^ fe . 
From Eqs. (17 19), we find the spectral decomposition of o~ z as 

^01+^02 . . / 7i + li\ , f woi - ^02 . if A 2 , ooA 

w ± = o ^ l — o — ) ± \ \ o h o(7i - 74) - 7273(22) 



where Re(u;±) are the main oscillation frequencies of the a z dynamics, and Im(u;±) gives 
their respective damping rates. 

Figure [2] presents a typical time evolution for a z for the resonant case using Eq. 
(22). Once again, by comparing it to the dynamics presented by a larger Hilbert space, 
iV = 18, we can verify a very good agreement of the results obtained using the scheme 
developed in section|2]. As one can observe, the TSS dynamics presents a rich structure, 
with the presence of two main oscillation frequencies, which are due to the TSS-HO 
resonance. These results corroborated those obtained in [T^l IT5l ITT]. 



4. Conclusions 



In conclusion, we presented and discussed a new approximation scheme to describe 
the dynamics of the spin-boson model with a structured environment. We studied the 
proposed problem by exploring its exact mapping onto that of a TSS coupled to a single 
HO of frequency ujq, which itself interacts with a bath of oscillators. We showed that, 
in order to find a low-dimensional effective Hamiltonian for the TSS-HO system, the 
representation of displaced HO states, DSR, provides an appropriate basis to truncate 
its Hilbert space. For that, we defined a conditional displacement operator, where 
the shifting of the HO system was dependent of the state of the TSS. In so doing 
we needed to introduce an ad hoc displacement parameter s which should be chosen 
in such a way that the best effective Hamiltonian could be obtained. Invoking the 
physics of an atom-cavity system, we chose the parameter s such that the effective 
Hamiltonian had a Jaynes-Cummings form. By comparing our numerical results with 
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Figure 2. a z (t) dissipative dynamics for the resonant case, A = u>q, computed for 
an Ohmic spectral density, Eq. pL with k = 0.02, uj c — 10A and temperature 
ksT — 0.1HA. We assumed a a z — coordinate TSS-HO coupling, with constant coupling 
given by g — 0.3A. Two main oscillation frequencies ui±, Eq. (22 1, are observed in 
the dynamics, reflecting the hybridization of the TSS-HO states. Solid and dotted 
lines represent, respectively, the N = 18 dimensional Hilbert space and our N = 4 
truncated space in the DSR. 



those considering a larger Hilbert space, we showed that, in fact, the derived effective 
Hamiltonian gives an excellent approximation of the original problem for a vast range of 
its physical parameters even when the TSS and the HO are moderately coupled close to 
the resonance. Actually, we showed that in order to achieve the same level of precision 
of the low-dimensional DSR model one needs to, at least, double the dimensionality 
of Hilbert space of the system within the simple truncation approximation. Finally, 
assuming the regime of weak HO-bath coupling and low temperatures, we were able to 
analytically determine the dissipative TSS dynamics using a 3LS Redfield theory. 
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